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[T I ' Abstract: 
OO . 

CN . We discuss the spectral analysis of quarkonium states in a hot medium of deconfined 
quarks and gluons, and we show that such an analysis provides a way to determine the 

>■ ' thermal properties of the quark-gluon plasma. 

in ■ 
^ , 

! 1. Introduction 

O 

O ' Since some thirty years it is known that besides the almost massless up and down quarks 

^ . of the everyday world, and the still relatively light strange quarks required to account for 
P^! the strange mesons and hyperons observed in hadron collisions, there are heavy quarks at 

the other end of the scale, whose bare masses alone are larger than those of most of the 
normal hadrons. These heavy quarks first showed up in the discovery of the J/ip meson 
[T], of mass of 3.1 GeV; it is a bound state of a charm quark (c) and its antiquark (c), 
each having a mass of some 1.2-1.5 GeV. On the next level there is the T meson ^j, with 
a mass of about 9.5 GeV, made up of a bottom or beauty quark-antiquark pair (bb), with 
5^ I each quark here having a mass around 4.5 - 4.8 GeV. Both charm and bottom quarks 
can of course also bind with normal light quarks, giving rise to open charm {D) and open 
beauty {B) mesons. The lightest of these 'light-heavy' mesons have masses of about 1.9 
GeV and 5.3 GeV, respectively. 

The bound states of a heavy quark Q and its antiquark Q are generally referred to as 
quarkonia. Besides the initially discovered vector ground states J/ip and T, both the cc 
and the bb systems give rise to a number of other stable bound states of different quantum 
numbers. They are stable in the sense that their mass is less than that of two light-heavy 
mesons, so that strong decays are forbidden. The measured stable charmonium spectrum 
contains the IS scalar rjc and vector J/ip, three IP states Xc (scalar, vector and tensor), 
and the 25* vector state ijj', whose mass is just below the open charm threshold. There 
are further charmonium states above the ip'] these can decay into DD pairs, and we shall 
here restrict our considerations only to quarkonia stable under strong interactions. 
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* Dedicated to Adriano Di Giacomo on the occasion of his 70th birthday 
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The study of quarkonia has played and continues to play a major role in many aspects 
of QCD: in spectroscopy, in production and decay, and last but not least as probe of 
hot QCD media And in the development of our understanding of and our feeling for 
the charm and the beauty of strong interactions, Adriano Di Giacomo has played and 
certainly will continue to play a major role. Neither my competence nor the space and 
time available here allow me to give an overview of these contributions; so let me just cite 
a few recent ones which are of particular importance to what I want to discuss here IHEj- 

Quarkonia are rather unusual hadrons. The masses of the light hadrons, in particular 
those of the non-strange mesons and baryons, arise almost entirely from the interaction 
energy of their nearly massless quark constituents. In contrast, the quarkonium masses are 
largely determined by the bare charm and bottom quark masses. These large quark masses 
allow a very straightforward calculation of many basic quarkonium properties, using non- 
relativistic potential theory. It is found that, in particular, the ground states and the 
lower excitation levels of quarkonia are very much smaller than the normal hadrons, and 
that they are very tightly bound. Now deconfinement is a matter of scales: when the 
separation between normal hadrons becomes much less than the size of these hadrons, 
they melt to form the quark-gluon plasma. What happens at this point to the much 
smaller quarkonia? When do they become dissociated? That is the main question we 
want to address here. We shall show that the disappearance of specific quarkonia signals 
the presence of a deconfined medium of a specific temperature [Hj. Thus the study of the 
quarkonium spectrum in a given medium is akin to the spectral analysis of stellar media, 
where the presence or absence of specific excitation lines allows a determination of the 
temperature of the stellar interior. 

We had defined quarkonia as bound states of heavy quarks which are stable under strong 
decay, i.e., Mcc < 2M£) for charmonia and M^^i < 2Mb for bottomonia. Since the quarks 
are heavy, with rric — 1.2 — 1.5 GeV for the charm and rrib ~ 4.5 — 4.8 GeV for the bottom 
quark, quarkonium spectroscopy can be studied quite well in non-relativistic potential 
theory. The simplest ("Cornell") confining potential [71 for a QQ at separation distance 
r has the form 

V(r)=ar~- (1) 
r 

with a string tension a ~ 0.2 GeV^ and a Coulomb-like term with a gauge coupling 
a ~ 7r/12. The corresponding Schrodinger equation 

/2m,-— V2 + \/(r)U,(r) = Mi<l>i(r) (2) 
I rric J 

then determines the bound state masses Mj and the wave functions $j(r), and with 

(r2) = J d'r r2|$,(r)|V J dh \Mr)\\ (3) 

the latter in turn provide the (squared) average bound state "radii" , here defined as the 
QQ separation for the state in question. 

The solution of eq. Q gives in fact a very good account of the full (spin- averaged) quarko- 
nium spectroscopy, as seen in Table 1 jH] . The line labelled AM shows the differences be- 
tween the experimental and the calculated values; they are in all cases less than 1 %. The 
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input parameters for these results were rric = 1.25 GeV, = 4.65 GeV, ^/o^ = 0.445 GeV, 
a = 7r/12. We see that in particular the J/ip and the lower-lying bottomonium states 
are very tightly bound {2M£)^b — ^ ^^"^ very small spatial size (rg ^ 2rh ~ 2 
fm). What happens to them in a hot and dense medium? 



state 


J/^ 


Xc 


^' 


T 


Xh 




X'b 




mass [GeV] 


3.10 


3.53 


3.68 


9.46 


9.99 


10.02 


10.26 


10.36 


[GeV] 


0.64 


0.20 


0.05 


1.10 


0.67 


0.54 


0.31 


0.20 


AM [GeV] 


0.02 


-0.03 


0.03 


0.06 


-0.06 


-0.06 


-0.08 


-0.07 


ro [fm] 


0.50 


0.72 


0.90 


0.28 


0.44 


0.56 


0.68 


0.78 



Table 1: Quarkonium Spectroscopy from Non-Relativistic Potential Theory [H] 



2. Interaction Range and Colour Screening 

Consider a colour-singlet bound state of a heavy quark Q and its antiquark Q, put into 
the medium in such a way that we can measure the energy of the system as function of 
the QQ separation r (see Fig. H}. The quarks are assumed to be heavy so that they are 
static and any energy changes indicate changes in the binding energy. We consider first 
the case of vanishing baryon density; at T = 0, the box is therefore empty. 
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Figure 1: String breaking for a QQ system 



In vacuum, i.e., at T = 0, the free energy of the QQ pair is assumed to have the string 
form [Zj 

F{r) ~ ar (4) 

where a ~ 0.16 GeV^ is the string tension as determined in the spectroscopy of heavy 
quark resonances (charmonium and bottomonium states). Thus F[r) increases with sep- 
aration distance; but when it reaches the value of a pair of dressed light quarks (about 
the mass of a p meson), it becomes energetically favorable to produce a qq pair from the 
vacuum, break the string and form two light-heavy mesons {Qq) and {Qq)- These can 
now be separated arbitrarily far without changing the energy of the system (Fig. Q). 
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The string breaking energy for charm quarks is found to be 

Fq = 2{Md - rUc) ~ 1.2 GeV; (5) 
for bottom quarks, one obtains the same value, 

Fq = 2{Mb - rrib) ~ 1.2 GeV, (6) 

using in both cases the quark mass values obtained in the solution leading to Table 3. 
Hence the onset of string breaking is evidently a property of the vacuum as a medium. It 
occurs when the two heavy quarks are separated by a distance 

ro ~ 1.2 GeV/a ~ 1.5 fm, (7) 

independent of the mass of the (heavy) quarks connected by the string. 

If we heat the system to get T > 0, the medium begins to contain light mesons, and 
the large distance QQ potential F{oo,T) decreases, since we can use these light hadrons 
to achieve an earlier string breaking through a kind of flip-flop recoupling of quark con- 
stituents jnj, resulting in an effective screening of the interquark force (see Fig. |2]). Near 
the deconfinement point, the hadron density increases rapidly, and hence the recoupling 
dissociation becomes much more effective, causing a considerable decrease of F{oo,T). 




Figure 2: In- medium string breaking through recoupling 

A further increase of T will eventually bring the medium to the deconfinement point 
Tc, at which chiral symmetry restoration causes a rather abrupt drop of the light quark 
dressing (equivalently, of the constituent quark mass), increasing strongly the density of 
constituents. As a consequence, F(oo,T) now continues to drop sharply. Above Tc, light 
quarks and gluons become deconfined colour charges, and this quark-gluon plasma leads to 
a colour screening, which limits the range of the strong interaction. The colour screening 
radius ro, which determines this range, is inversely proportional to the density of charges, 
so that it decreases with increasing temperature. As a result, the QQ interaction becomes 
more and more short-ranged. 

In summary, starting from T = 0, the QQ probe first tests vacuum string breaking, then 
a screening-like dissociation through recoupling of constituent quarks, and finally genuine 
colour screening. In Fig. El we show the behaviour obtained in full two-flavour QCD for 
the colour-singlet QQ free energy as a function of r for different T [TO] . 

It is evident in Fig. |21 that the asymptotic value F{oo,T), i.e., the energy needed to 
separate the QQ pair, decreases with increasing temperature, as does the separation 
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Figure 3: The colour singlet QQ free energy F{r, T) vs. r at different T JU] 

distance at which "the string breaks" . For the moment we consider the latter to be defined 
by the point beyond which the free energy remains constant within errors, returning in 
section 4.3 to a more precise definition. The behaviour of both quantities is shown in Fig. 
0] Deconfinement is thus reflected very clearly in the temperature behaviour of the heavy 
quark potential: both the string breaking energy and the interaction range drop sharply 
around Tc. The latter decreases from hadronic size in the confinement region to much 
smaller values in the deconfined medium, where colour screening is operative. 



■o<T) [I'mJ 




0.5 1.0 1.5 T/Tc 0.5 1.0 1.5 

Figure 4: String breaking potential and interaction range at different temperatures 

The in-medium behaviour of heavy quark bound states thus does serve quite well as probe 
of the state of matter in QCD thermodynamics. We had so far just considered QQ bound 
states in general. Let us now turn to a specific state such as the J/ip- What happens when 
the range of the binding force becomes smaller than the radius of the state? Since the c 
and the c can now no longer see each other, the J/tp must dissociate for temperatures 
above this point. Hence the dissociation points of the different quarkonium states provide 
a way to measure the temperature of the medium. The effect is illustrated schematically 
in Fig. 13 showing how with increasing temperature the different charmonium states 
"melt" sequentially as function of their binding strength; the most loosely bound state 
disappears first, the ground state last. 

Moreover, since finite temperature lattice QCD also provides the temperature dependence 
of the energy density the melting of the different charmonia or bottomonia can be 
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specified as well in terms of e. In FiglHl we illustrate this, combining e(T) with the force 
radii shown in Fig. EJ It is evident that although ip' and Xc are both expected to melt 
around Tc, the corresponding dissociation energy densities could still differ. 
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Figure 5: Charmonium spectra at different temperatures 




Figure 6: Charmonium dissociation vs. temperature and energy density 



To make these considerations quantitative, we thus have to find a way to determine the 
in-medium melting points of the different quarkonium states. This problem has been 
addressed in three different approaches: 

• Model the heavy quark potential as function of the temperature, V{r, T), and solve 
the resulting Schrodinger equation 

• Determine the internal energy [/(r, T) of a QQ pair at separation distance r from 
lattice results for the corresponding free energy F{r,T), using the thermodynamic 
relation 

U^r, T).-T^ i^-^l^) ^ Fir, T) - T i?^) . (8) 

and solve the Schrodinger equation with V{r, T) = U (r, T) as the binding potential. 

• Calculate the quarkonium spectrum directly in finite temperature lattice QCD. 
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Clearly the last is the only model-independent way, and it will in the long run provide 
the decisive determination. However, the direct lattice study of charmonium spectra has 
become possible only quite recently, and so far most results are obtained in quenched QCD 
(no dynamical light quarks); corresponding studies for bottomonia are still more difficult. 
Hence much of what is known so far is based on Schrodinger equation studies with different 
model inputs. To illustrate the model-dependence of the dissociation parameters, we first 
cite some early work using different models for the temperature dependence of V{r,T), 
then some recent studies based on lattice results for F{r,T), and finally summarize the 
present state of direct lattice calculations of charmonia in finite temperature media. 

3. Potential Model Studies 

The first quantitative studies of finite temperature charmonium dissociation [12j were 
based on screening in the form obtained in one-dimensional QED, the so-called Schwinger 
model. The confining part of the Cornell potential (0), V{r) ~ or, is the solution of the 
Laplace equation in one space dimension. In this case, Debye-screening leads to ^3] 



where /i(T) denotes the screening mass (inverse Debye radius) for the medium at tem- 
perature T. This form reproduces at least qualitatively the convergence to a finite large 
distance value V{oo,T) = cr//i(T), and since /i(T) increases with T, it also gives the 
expected decrease of the potential with increasing temperature. Combining this with the 
usual Debye screening for the 1/r part of eq. (P) then leads to 



for the screened Cornell potential. In |12J, the screening mass was assumed to have the 
form fi{T) ~ 4 T, as obtained in first lattice estimates of screening in high temperature 
SU (N) gauge theory. Solving the Schrodinger equation with these inputs, one found that 
both the ip' and the Xc become dissociated essentially at T ~ Tc, while the J/ip persisted 
up to about 1.2 Tc. Note that as function of the energy density e ~ T^, this meant that 
the J/ijj really survives up to much higher e. 

This approach has two basic shortcomings: 

• The Schwinger form Q corresponds to the screening of ar in one space dimension; 
the correct result in three space dimensions is different 

• The screening mass /i(r) is assumed in its high energy form; lattice studies show 
today that its behaviour near is quite different [T^ . 

While the overall behaviour of this approach provides some first insight into the problem, 
quantitative aspects require a more careful treatment. 

When lattice results for the heavy quark free energy as function of the temperature first 
became available, an alternative description appeared jTH]- It assumed that in the ther- 
modynamic relation (jH)) the entropy term —T{dF/dT) could be neglected, thus equating 
binding potential and free energy. 




(9) 




(10) 



V{r,T) = F{r,T) 



T{dF/dT) ^ F(r,T). 



(11) 
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Using this potential in the Schrodinger equation (0) specifies the temperature dependence 
of the different charmonium masses. On the other hand, the large distance limit of V{r, T) 
determines the temperature variation of the open charm meson D, 

2Mb (T) ~ 2me + 1/(00, T) (12) 

In fig. we compare the resulting open and hidden charm masses as function of tem- 
perature. It is seen that the mass falls below 2 Md around 0.2 Tc, that of the Xc at 
about 0.8 Tf.; hence these states disappear by strong decay at the quoted temperatures. 
Only the ground state J/ip survives up to Tc and perhaps slightly above; the lattice data 
available at the time did not extend above Tc, so further predictions were not possible. 

3.8 



^ 3.4 
> 

O 

3.0 
2.6 

"V- 0.5 Xc 1 

Figure 7: Temperature dependence of open and hidden charm masses jTHI 




The main shortcoming of this approach is quite evident. The neglect of the entropy term 
in the potential reduces V{r,T) and hence the binding. As a result, the D mass drops 
faster with temperature than that of the charmonium states, and it is this effect which 
leads to the early charmonium dissociation. Moreover, in the lattice studies used here, 
only the colour averaged free energy was calculated, which leads to a further reduction of 
the binding force. 

We conclude from these attempts that for a quantitative potential theory study, the 
free energy has to be formulated in the correct three-dimensional screened Cornell form, 
and it then has to be checked against the space- and temperature-dependence of the 
corresponding colour singlet quantity obtained in lattice QCD. 

4. Screening Theory 

The modification of the interaction between two charges immersed in a dilute medium of 
charged constituents is provided by Debye-Hiickel theory, which for the Coulomb potential 
in three space dimensions leads to the well-known Debye screening, 

I ^ y^^^ (13) 

where rij = defines the screening radius [in]. Screening can be evaluated more 
generally ^3] for a given free energy F{r) ~ r"^ in c? space dimensions, with an arbitrary 
number q. We shall here apply this to the two terms of the Cornell form, with q = 1 for 
the string term, q = —1 for the gauge term, in d = 3 space dimensions |14j. 
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We thus assume that the screening effect can be calculated separately for each term, so 
that the screened free energy becomes 



a 



F{r,T) = F,(r,T) + F,(r,T) = ar/,(r,T) - -/e(r,T). 

r 

The screening functions fs{r,T) and fc{r,T) must satisfy 

/,(r,T) = Ur,T) = l for T ^ 0, 
/s(r,T) = Ur,T) = l forr^O, 

since at T = there is no medium, while in the short- distance limit T^^ ^ r 



a 
r 



fir 



medium has no effect. The resulting forms are 

Fcir,T) 

for the gauge term, and 

F.(rX) = -^ r(l/4) 



23/2r(3/4) 23/41(3/4) 



(14) 



(15) 
0, the 

(16) 
(17) 



for the string term. The first term in eq. (jl7|) gives the constant large distance limit due 
to colour screening; the second provides a Gaussian cut-off in x = fir, since Ki/4{x'^) ~ 
exp{— x^)); this is in contrast to the exponential cut-off given by the Schwinger form ()10|). 
but in accord with string breaking arguments of Di Giacomo et al. 





(a) (b) 
Figure 8: Screening fits to the QQ free energy F{r, T) for T >Tc (left) and T < (right) 



At temperatures T > Tc, when the medium really consists of unbound colour charges, we 
thus expect the free energy F{r,T) to have the form ()14|) . with the two screened terms 
given by eqs. (fT^ and (fTTj) . In Fig. |S^, it is seen that the results for the colour singlet 
free energy calculated in two-flavour QCD ^7] are indeed described very well by this 
form, with rric = 1.25 GeV and ^/a = 0.445 GeV, as above. The only parameter to be 
determined, the screening mass fi, is shown in Fig. El and as expected, it flrst increases 
rapidly in the transition region and then turns into the perturbative form /i ~ T. 
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The behaviour of QQ binding in a plasma of unbound quarks and gluons is thus well 
described by colour screening. Such a description is in fact found to work well also for 
T < Tc, when quark recombination leads to an effective screening-like reduction of the 
interaction range, provided one allows higher order contributions in x = /xr in the Bessel 
function i^i/4(a;^) governing string screening ^3]. The resulting fit to the two-flavour 
colour singlet free energy below is shown in Fig. |S|d, using Kin^x^ + {l/2)x^) in eq. 
fll7j) : the corresponding values for the screening mass are included in Fig. IHl 

1.6 
1.4 
1.2 

1 

0.8 
0.6 
0.4 
0.2 

°0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

Figure 9: The screening mass /i(T) vs. T [EllHl 

With the free energy F(r, T) of a heavy quark- antiquark pair given in terms of the screen- 
ing form obtained from Debye-Hiickel theory, the internal energy f/(r, T) can now be 
obtained through the thermodynamic relation (jH}, and this then provides the binding 
potential V{r,T) for the temperature dependent version of the Schrodinger equation 
The resulting solution then specifies the temperature-dependence of charmonium binding 
as based on the correct heavy quark potential |H]. Let us briefly comment on the ther- 
modynamic basis of this approach. The pressure P of a thermodynamic system is given 
by the free energy, P = —F = —U + TS; it is determined by the kinetic energy TS at 
temperature T and entropy S(T), reduced by the potential energy U(T) between the con- 
stituents. In our case, all quantities give the difference between a thermodynamic system 
containing a QQ pair and the corresponding system without such a pair. The potential 
energy of the QQ pair, due both to the attraction of Q and Q and to the modification 
which the pair causes to the internal energy of the other constituents of the medium, is 
therefore given by U. To determine the dissociation points for the different quarkonium 
states, we thus have to solve the Schrodinger equation ^ with V{r,T) = U{r,T). 

From eqs. (fTH|) and (fT7|) we get for the QQ potential 

V{r,T) = V{oo,T) + V{r,T), (18) 

with 

V{oo,T)=ci an + T—[ci— + a], (19) 

/i al /i^ 

and Ci = r(l/4)/2^/^r(3/4), and where K(r, T) contains the part of the potential which 
vanishes for r oo. The behaviour of V{oo,T) as function of the temperature is shown 
in Fig. Uni It measures (twice) the energy of the cloud of light quarks and gluons around 
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an isolated heavy quark, of an extension determined by the screening radius, relative to 
the energy contained in such a cloud of the same size without a heavy quark. This energy 
difference arises from the interaction of the heavy quark with light quarks and gluons of 
the medium, and from the modification of the interaction between the light constituents 
themselves, caused by the presence of the heavy charge. - The behaviour of V{r, T) 
is shown in Fig. ^2 for three different values of the temperature. It is seen that with 
increasing T, screening reduces the range of the potential. 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
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Figure 10: The large distance limit of 
the quarkonium potential V{r,T) 



Figure 11: Variation of V(r, T) with r 
for different T |Hj 



The relevant Schrodinger equation now becomes 

/±V2-\/(r,T)U,(r) = AE,(T)$,(r) (20) 

where 

AEi{T) = V{oo, T)- Mi- 2mc (21) 

is the binding energy of charmonium state i at temperature T. When it vanishes, the 
bound state i no longer exists, so that AEi{T) = determines the dissociation temper- 
ature Tj for that state. The temperature enters only through the T-dependence of the 
screening mass /i(T), as obtained from the analysis of the lattice results for F{r,T). - 
Equivalently, the divergence of the binding radii, rj(T) oo, can be used to define the 
dissociation points T,. 

The same formalism, with = 4.65 GeV replacing rric, leads to the bottomonium dis- 
sociation points. The combined quarkonium results are listed in Table 2. They agree 
quite well with those obtained in a very similar study based on corresponding free en- 
ergies obtained in quenched lattice QCD [TH], indicating that above Tc gluonic effects 
dominate. Using a parametrically generalized screened Coulomb potential obtained from 
lattice QCD results also leads to very compatible results for the Nf = 2 and quenched 
cases P^. We recall here that the main underlying change, which is responsible for the 
much higher dissociation temperatures for the quarkonium ground states, is the use of 
the full internal energy (jH}, including the entropy term, as potential in the Schrodinger 
equation: this makes the binding much stronger. 
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We should note, however, that in all such potential studies it is not so clear what binding 
energies of less than a few MeV or bound state radii of several fermi can mean in a medium 
whose temperature is above 200 MeV and which leads to screening radii of less than 0.5 
fm. In such a situation, thermal activation [201 can easily dissciate the bound state. 



state 


J/^(15) 


Xc(lP) 


^'(25) 


T(15) 


Xb{lP) 


T{2S) 


Xb{2P) 


T(35) 




2.10 


1.16 


1.12 


> 4.0 


1.76 


1.60 


1.19 


1.17 



Table 2: Quarkonium Dissociation Temperatures [Hj 
5. Charmonium Correlators 

The direct spectral analysis of charmonia in finite temperature lattice has come within 
reach only in very recent years [21] ■ It is possible now to evaluate the correlation functions 
Gh{t, T) for hadronic quantum number channels H in terms of the Euclidean time r and 
the temperature T. These correlation functions are directly related to the corresponding 
spectral function aH{M,T), which describe the distribution in mass M at temperature 
T for the channel in question. In Fig. ^[ schematic results at different temperatures are 
shown for the J/ ip and the Xc channels. It is seen that the spectrum for the ground state 
J /ifj remains essentially unchanged even at 1.5 Tc. At 3 Tc, however, it has disappeared; 
the remaining spectrum is that of the cc continuum of J /if) quantum numbers at that 
temperature. In contrast, the Xc is already absent at 1.1 Tc, with only the corresponding 
continuum present. 



J/v Zc 




2 3 4 2 3 

M TGeVl M TGeVl 



Figure 12: J/ifj and Xc spectral functions at different temperatures 

These results are clearly very promising: they show that in a foreseeable future, the 
dissociation parameters of quarkonia can be determined ah initio in lattice QCD. For 
the moment, however, they remain indicative only, since the underlying calculations were 
generally performed in quenched QCD, i.e., without dynamical quark loops. Since such 
loops are crucial in the break-up of quarkonia into light-heavy mesons, final results require 
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calculations in full QCD. Some first calculations in two-flavour QCD have just appeared 
[22] and support the late dissociation of the J/ip. The widths of the observed spectral 
signals are at present determined by the precision of the lattice calculations; to study 
the actual physical widths, much higher precision is needed. Finally, one has so far only 
first signals at a few selected points; a temperature scan also requires higher performance 
computational facilities. Since the next generation of computers, in the multi-Terafiops 
range, is presently going into operation, the next years should bring the desired results. 
So far, in view of the mentioned uncertainties in both approaches, the results from direct 
lattice studies and those from the potential model calculations of the previous section 
appear quite compatible. 

6. Conclusions 

The theoretical analysis of the in-medium behaviour of quarkonia has greatly advanced in 
the past decade. Potential model studies based on lattice results for the colour-singlet free 
energy appear to be converging to results from direct lattice calculations, and within a 
few more years, the dissociation temperatures for the different quarkonium states should 
be known precisely. Through corresponding calculations of the QCD equation of state, 
these temperatures provide the energy density values at which the dissociation occurs. In 
statistical QCD, quarkonia thus allow a spectral analysis of the quark-gluon plasma, as 
many of the observed features as possible. 

Finally, it seems worthwhile to note that experimental measurements of the relative dis- 
sociation points of the different quarkonium states could be possible in the study of high 
energy nucleus-nucleus collisions. If they succeed, this might in fact be a unique chance 
to compare quantitative ah initio QCD predictions directly to data. 
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